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• We study the phase-ordering kinetics of the one-dimensional Heisenberg model with conserved 

' order parameter, by means of scaling arguments and numerical simulations. We find a rich dynamical 

pattern with a regime characterized by two distinct growing lengths. Spins are found to be coplanar 
over regions of a typical size Lv{t), while inside these regions smooth rotations associated to a smaller 
' length Lc{t) are observed. Two different and coexisting ordering mechanisms are associated to these 

^ , lengths, leading to different growth laws Lv{t) ~ t^'''^ and Lc{t) ~ t^'''^ violating dynamical scaling. 
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I. INTRODUCTION 



After quenching a ferromagnetic system to a low temperature phase, relaxation towards the new equilibrium state is 
realized by a progressive phase-ordering |T] . The specific mechanisms involved in the coarsening phenomenon depend 
on the presence and the nature of topological defects. In d-dimensional systems described by and 0{N) vector order 
parameter, topological defects are unstable for iV > d -I- 1 (for N — d + 1 peculiar defects as textures [2| may be 
Q ] present). Therefore, in the asymptotic regime when all defects have disappeared, the dynamics is solely driven by 
'— the reduction of the excess energy related to the smooth rotations of the order parameter. In contrast, systems 
^ ' with N < d are characterized by the presence of stable defects whose presence influences the dynamics in the whole 
^ , phase-ordering stage. In particular, when N = d defects are localized and ordering occurs by mutual defect-antidefect 
T^-j- ■ annihilation. This is the case of the Ising chain, where up and down domains are separated by point-like interfaces 
performing random walks. 

Generally, the late stage is characterized by dynamical scaling [l|, 0] . This implies that a single characteristic length 
L{t) can be associated to the development of order in such a way that configurations of the system are statistically 
independent of time when lengths are measured in units of L{t). The characteristic length usually has apower law 
^ growth L{t) oc t^/^ . In systems with a conserved order parameter (COP) one generally finds z = 3[l, [H|orz = 4 
On . [H, Q for N — 1 and > 1 respectively. 

' For systems at or below the lower critical dimension d^, such as the Ising chain, a true asymptotic phase-ordering 
^ , can only be observed in quenches to T = 0. However, if quenches to a relatively low temperature are performed, one 
observes an initial transient regime (but very long lasting when T is small) where the dynamics is indistinguishable 
from that at T=0. This regime lasts until L{t) has grown comparable to the equilibrium coherence length £,{T). 
In this Article, we investigate the phase-ordering kinetics of the one-dimensional Heisenberg model {N — 3) with 
■ ■ ' COP quenched to a low temperature T. We show that the dynamics is much richer than what one would naively 
expect. This is due to the formation, in an early stage, of couples of parallel spins, separating regions in which the spins 
are coplanar. These parallel spins act as point-like defects, as it will be explained in Sec. IIIII Their presence provides 
an analogy between the Heisenberg and the Ising chain, where couples of parallel spins and regions of coplanarity in 
the former model correspond to interfaces and domains in the latter. The analogy is not only formal, but is reflected 
in the kinetics: in a first stage (whose duration however diverges in the T — s- limit) regions of coplanar spins coarsen 
similarly to the domains of the Ising model, their typical length growing as 

Lv{t)^t^l\ (1) 



In this regime the number of defects is reduced only by a mechanism which recalls the annihilation of the interfaces in 
the Ising model. Simultaneously, smooth rotations of the spins, typical of vectorial systems, occur inside the regions 
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of coplanarity. The coherence of the spins inside these regions extends over a length Lc{t) ^ Lv{t) increasing as 

Lc{t)^t''\ (2) 

The existence of two growing length, associated to different ordering mechanisms, produces the breakdown of dynam- 
ical scaling. In the analogy between the Heisenberg and the Ising chain, a notable difference must be stressed. While 
in the latter interfaces are stable defects which can only be removed by mutual annihilation, defects in the former 
are unstable. Namely, after a first stage of coarsening of coplanar regions, the defects spontaneously decay due to 
thermal fluctuations. The typical lifetime of the defects being limited by temperature, coarsening of the coplanar 
regions persists up to very long times in deep quenches. After that, defects disappear and the system finally enters 
a late stage where smooth spin rotations remains the only mechanism at work, until equilibration is attained when 
Lc{t) ^ £.{T)- In this regime dynamical scaling is restored with exponent z = 4, as expected for a system with 
N > d + I. In the small temperature limit the duration of the two regimes, with and without dynamical scaling, are 
comparable. 

This paper is organized as follows: In Sec.|lT]we introduce the model and define the observable quantities that will 
be considered. In Sec. IIIII we describe the main features of the dynamics in the different regimes, compute the value 
of the exponents and of other quantities by means of scaling arguments, and compare our results with the outcome 
of numerical simulations of the model. A summary and the conclusions are contained in Sec. IIVI 



II. MODEL AND OBSERVABLES 



The Heisenberg model is defined by the Hamiltonian 

H[(j] = = - J • - 1) = - J ^(cos(t>i - 1), (3) 

where is the local energy density, Ui is a three-components unit vector spin, i — 1, A/" are the sites on a one- 
dimensional lattice and (j)i is the angle between di and Ui+i. We will assume J = 1 and the Boltzmann constant 
kB = l. 

The equilibrium properties of the model are exactly known 7]. This system is ergodic except at T = 0. At any 
finite temperature the state is disordered with a vanishing magnetization and internal energy (per spin) Eeq{T) = 
T — coth(l/r) + 1 with the low temperature expansion Ef.q{T) T. The correlation function Ceq{r) = {uiCTi^r) = 
[1 — Ef.q{T)]'' decays exponentially over a coherence length ^(T) that diverges in the T ^ limit. 

Concerning dynamics, at equilibrium the model possesses an intrinsic kinetics where the energy and the magneti- 
zation are conserved Q . Studies Q of this intrinsic dynamics have evidenced the failure of the equilibrium scaling 
symmetry close to T = 0. Experimental work supports this picture [lo| . 

In this paper we consider a different situation, where the system is quenched from an high temperature configuration 
to a low temperature T, transferring the excess energy to a heat bath coupled to it. In order to model this, we consider 
a dynamics where two neighboring spins of a configuration [a] are randomly chosen at each time-step and 

then they are updated to provided the local magnetization is conserved, namely Si = Oi + (t^+i = a[ + ct^^i. 

Notice that magnetization is conserved at the local level with this rule. Due to the conservation law, the spins (7^, iJi+i 
can only rotate rigidly around their sum s^, as shown in Fig. [TJ We consider the heat-bath transition rates [a] — *■ [(?'] 
satisfying detailed balance 

wAa']^W-'e^^[-^y (4) 

where Wi = J da[da[^i5{ai + Si+i ~ s*;) exp(— _ff [(t']/T). Heat-bath transition rates provide a particularly fast and 

efficient dynamics with respect to other (i.e. Metropolis) choices [ll|. Let us denote di = Oi+i — Oi, d[ — — a[ 
and Pi the projection of (7^+2 ^ o^i-i on the plane H perpendicular to Si (see Fig. [T]). Any move, involving the couple 
(Ji and (7i4-i, can be described by a rotation in the plane H from di to d'j. In this framework the angles between pi 
and di (d'i), denoted as 6i {0[), fully parametrize the dynamics and the transition rate ^ can be rewritten as 

/Q'\ W-l f d^p^{cOSe'i- 1)\ 

expl — 1, (5) 

where di = \di\ = \di\, Pi = \pi\. Von Neuman rejection method [T2j allows to efficiently generate 9'^ according to the 
transition rates ^ . Notice that in a move the typical deviations from the lowest energy configuration {6'^ = 0) are of 
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FIG. 1: Schematic representation of the four spins tji, ai+i, 5^+2. 



order 

2T 

cos^^-l^— . (6) 
diPi 

We consider a system initially prepared in an high temperature uncorrelated state, with Yl^=i — ^i ^"^^ then 
quenched, at time i = 0, to a lower final temperature T . As already mentioned in the Introduction, the dynamics of 
systems at or below the critical dimension, such as the one considered here, is characterized by an initial transient 
where the system orders [l3| as in a quench to T = 0. The characteristic length of ordered regions grows in time until, 
at time Teq{T) it becomes comparable to £,{T). At this point the final equilibrium state at T is entered and phase- 
ordering ends. If the system is quenched to a sufhciently low temperature, since £,{T) is very large, the phase-ordering 
kinetics extends over a huge time window t < Teq{T). 

Characteristic lengths, and scaling properties, can be studied from the knowledge of the two-points equal time 
correlation function 

C{r,t) = {a,{t)-a,+r{t)), (7) 

where (...) means an ensemble average, namely taken over different initial conditions and thermal histories. Due to 
space homogeneity, C{r,t) does not depend on i. Dynamical scaling would imply 

C(r,i)=c(x), (8) 

where x = r/L{t). For systems with with an 0{N) vector order parameter the Bray-Puri-Toyoki TF] behavior 

1 - c(x) ~ x^ (9) 



4 



is observed for small x. In the scalar case, this behavior reduces to the Porods law [l6| 

l~c{x)r^X, (10) 

which is generally expected in systems with sharp interfaces. From Eq. ([5]) one can extract a quantity Lc{t) propor- 
tional to L{t) from the condition 

C[Lcit),t]^^, (11) 

namely as the half- height width of C{r, t). In the following we will also consider the correlation 

V{r,t) = {v,it)-v,+r{t)). (12) 

where the unit vectors 



(13) 



identify the planes formed by neighboring spins, hence V{r.t) represents the correlations between these planes. When 
scaling holds, V{r,t) behaves similarly to C{r,t), namely 

Vir,t) = v{x), (14) 

and, defining Lv{t) analogously to Lc{t) in Eq.(' [TT|) . also Lv{t) oc L{t). 



III. THE DYNAMICS: SCALING ARGUMENTS AND NUMERICAL RESULTS 



In the following we will discuss the main features of the dynamical process, by means of scaling arguments and 
numerical simulations. These are performed on a string of 8000 spins with periodic boundary conditions cr^+i — a^. 
We have checked that with this size our simulations are free from finite size effects. An average over 5 — 10 realizations 
is made for each simulation. In the limit of low temperatures, to which we are interested in, the kinetics becomes 
very slow since, as it will be discussed below, time rescales as T"^/^; simulations are therefore quite time consuming 
in this region. 

The dynamics of a low temperature quench is characterized by different subsequent regimes, which are separately 
discussed below. 



A. Pinning (Quenches to T = 0). 



The kinetics of the quench to T = is determined by the existence of frozen states where the system gets trapped 
after a while. The nature of these states can be understood by looking at the first line (denoted as time i) of Figure [2l 
The two spins on sites i,i -I- 1 are parallel, hence no move involving this couple can be done, since the angle di is 
not defined. On the left and on the right of the aligned spins there are regions where the spins are coplanar, and 
hence the Vi are parallel and point along certain directions, denoted by a, etc., which can be considered as different 
phases of the system. These phases extend up to another couple of parallel spins (not shown in the figure). As it will 
be shown in Sec lIII Cl in quenches to finite temperatures the system depins after a while and these coplanar regions 
coarsen much in the same way as equilibrium phases do in usual coarsening systems. Due to this analogy, the term 
phases is used also here. However it must be noticed that regions of coplanarity are not equilibrium phases, because 
in equilibrium spins are coplanar AND aligned, while here they typically rotate (see Fig. [21), as will be discussed 
in Sec. IIII CI For the following discussions, we also introduce the terminology of distance between two phases a, /3, 
related to the angle by which the vectors Vi of the phase a must rotate in order to align with those of the phase f3. 
In this sense we will also talk of nearby or distant phases. 

As anticipated in Sec. [H the impossibility to eliminate couples of parallel spins by means of local rearrangements, 
involving tTi, (7^+1 alone, make them reminiscent of topological defects. Actually these spins represent real topological 
defects, as is readily seen by considering the representation in terms of Vi instead of (7^. With this description, phases 
are domains in a strict sense, namely Vi is constant in the interior of coplanar regions. They will be denoted as 
domains in the following, or planes, in view of the coplanarity of their spins, without further specification. Different 
domains are separated by sharp boundaries and, right on top of them, there is a localized defect where is not 



5 





FIG. 2: Schematic representation of spin configurations at three subsequent times t,t + l,t + 2. For each time the upper and 
lower lines represent the a and v configuration. Given a vector a (or v) with components (a, b, c), the arrow in the figure is a 
vector of components (a,fe, 0), namely the projection of a on the two-dimensional plane of the figure. The third component c 
can be read off by the constraint of unitary length of a, with the help of the unitary circles represented around each vector. 
The origin of a vector is marked with an heavy dot when c < (vector pointing behind the figure) , or with an open circle when 
c > 0. The meaning of the phases a, /3, 7 is discussed in the text. 



defined. Then, once the proper representation is considered, parallel spins qualify as defects in the usual sense. Since 
spins are coplanar inside the domains and parallel on a defect, it is clear that any move involving any couple of spins is 
forbidden at T = 0, and the dynamics is frozen on states like those depicted on the first line of Fig. [2l By identi fyin g 
aligned spins as defects separating domains, an analogy with the COP Ising model, which freezes as well at T = 17|, 
can be drawn. Let us mention, however, at least two main differences. First, the constraint on the motion of parallel 
spins in the Heisenberg model is related to the kinetic rule and not to T = 0. Second, as will be discussed in Sec. lIII Dl 
the defects in the Heisenberg model are unstable, although their lifetime diverges in the T — » limit. 

When a quench to T = is performed, the system starts reducing its energy by ordering the spins until some 
couples happen to be nearly parallel. In the meanwhile the n spins between two defects adapt themselves on a plane. 
Since n is a finite number this process can be accomplished in a finite time. At this point, the model gets trapped in 
one of the absorbing states discussed above. Notice that, in a situation as the one discussed here, Lv{t) and Lc{t) 
describe, respectively, the length of the domains and the coherence length of the spins in the bulk of the planes. The 
evolution of the model toward the pinned state can be studied by following the evolution of these lengths in the insets 
of Figs. [21 m In a quench to T = both these quantities initially grow but then saturate to a constant value when 
the system freezes. 
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FIG. 3: The quantity Lc{t) is plotted against rescaled time tT^^^ for different temperatures. The dot-dashed line is the 
expected behavior Lc{t) oc (tT^^^y^'^ after the depinning (see text). In the inset the same quantity is plotted against time in 
the case T = 0. 



B. Depinning 

When T 7^ but sufficiently small, the dynamics leading the system to the frozen state proceeds practically as in 
the case T = described in Sec. If 11 Al since it is entirely dominated by the moves which lower the energy which are 
not affected much by a small T. When the system is trapped in the absorbing state, however, activated moves can 
occur inside the domains if T > 0. According to Eq. since the system is still very disordered in this stage, both 
di and pi are on average large, and the typical value 9 of the angle 9i (we will use the over-bar in the following to 
denote the typical value of a quantity) is small. Notice, in fact, that even for the smallest temperatures considered 
in the simulations reported in Figs. ISBl pinning is never complete: Lc(t) and Lv{t) keep slowly increasing because 
couples of spins are tiny twisted. This mechanism eventually depins the system from the absorbing state, restoring 
the dynamics, as shown in Figs. 13141 after a characteristic time Tp{T). From Eq. ([6]), assuming 9 to be small, one has 

cos9 — 1 ~ ~ 2T/{dp) and hence 9 T^/^. Since 9 is tiny, a number n{T) n/9 oi these moves is needed in order 
to produce an appreciable decorrelation (of an angle of order, say, tt) with respect to the pinned state and to restart 
the dynamics. Therefore we find 

Tp(T) ~ (15) 

According to this result, for sufficiently low temperatures the curves for Lc{t) or iy(i) (and, more generally of any 
time dependent observable), should collapse after the pinning stage when plotted against tT^/'^. Figs. [31 [5] show 
that this is indeed quite well verified over 10 decades in temperature (with some deviations that will be discussed in 

Sec, mm . 
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FIG. 4: The quantity Lv{t) is plotted against rescaled time tT^^'^ for different temperatures. The dot-dashed line is the 
expected behavior Lv{t) oc (tT^^'^Y^^ in the first phase-ordering regime (see text). In the inset the same quantity is plotted 
against time in the case T = 0. 

C. First phase-ordering regime: presence of domains 

Here we give a schematic description of the microscopic kinetics in a first dynamical regime occurring after the 
depinning. In this regime spins evolve in such a way that the coplanarity of the spins inside the phases is preserved. 
Actually, the domains compete among themselves and grow much in the same way as the equilibrium phases in usual 
coarsening systems, as testified bye the increase of Lv{t) (see Fig. 14]). This regime lasts until the phenomenon of the 
breakdown of the plains, discussed in Sec. IIIID| occurs (when this happens Lv{t) stops growing and goes abruptly 
to zero. The end of this first dynamical stage can then be easily recognized by inspection of Fig. 2]). 

As we will explain below, in this regime two mechanisms are at work: The former is responsible for the coarsening 
of the domains, the latter is the phase-ordering of the spins inside the domains. Since these mechanisms are associated 
to two different growing lengths, dynamical scaling is not obeyed, as we will show explicitly. 

1. Coarsening of the domains 

In the limit of small T, among the moves discussed in Sec. IIII Bl those which produce the smallest energy increase 
are overwhelmingly favored. These are generally the moves involving the spins near the boundary of a plane, for 
instance those on sites i — 1, i in Fig. [21 The effect of this move occurring in the /3-phase, is the nucleation of a third 
phase, denoted by 7, as it is clear considering the direction of the vector Vi-i in the second line (denoted as time t + 1) 
of the figure (a spurious phase on site i — 2 is also generated, whose presence is, however, irrelevant). At this point 
the dynamics may proceed by rotating the spins (7i_2,i?i-i in order to make (7i_2, (Xi-i, and (li coplanar (third line, 
denoted as time t + 2, in the figure). In this way the new phase 7 may spread replacing the pre-existing phase (3. After 
the complete replacement of the phase [3 with 7 (when this occurs), if the energy of the system is increased (because 7 
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is more distant to the neighboring phase a (or the one on the left, not shown in the figure) than the original /3-phase, 
the 7 phase is quickly re-adsorbed by reversing the process. Conversely, if the energy of the system is decreased a 
new blocked state is reached, characterized by domains of more nearby phases. This mechanism provides a direction 
to the process, favoring on the average the formation of new phases for which a diminishing of the system energy 
occurs. At this point the process can start again with the activated nucleation of a new phase replacing 7 and so on 
repeatedly until one of the two adjacent phases (say a) extends over the original domain of the /? phase, increasing 
the typical size of the domains. 

The basic steps of this ordering process may recall what happens in the COP Ising model [13] . Actually, in both 
cases there are domains of different phases (di = ±1 in Ising, Vi = a, (3 etc. in the Heisenberg model) separated 
by sharp interfaces (see also the discussion on Porod tails below in this section). In both cases, the first step is 
the nucleation of a germ of another phase (the evaporation of a monomer in Ising) inside a domain of a preexistent 
phase. After nucleation, the kinetics proceeds by a random motion of the nucleated phase. This analogy is not only 
qualitative, since, as we show below, the growth exponent of Lv(t) is the same {z = 3). 

Actually this can be inferred by the following argument: After the first move, the position of the boundary between 
the new nucleated phase 7 and the remaining of the old f3 phase performs a random walk. Considering the long time 
regime, where Lv{t) is large, most of the times the boundary returns to its original position j. In this case the 7 
phase is re-adsorbed after the duration 5 ^ Ly{t) of the random walk. The probability that the interface moves a 
distance Lv{t) and hence the /3-phase is eliminated is proportional to 1/ Lv{t) [13]; for this reason the whole process 
must be repeated Ly(t) times in order to substitute the old phase (3. This requires a time 

T^(^Lv{t)5ozLv{tf. (16) 

This conjecture has been tested by means of numerical simulations, mimicking the evolution of a single domain as 
follows: We have prepared a domain of Ly initially coplanar spins, with a uniform rotation such that = e, Vz in the 
bulk of the domain. The boundary condition is made of two spins (on each side) lying on a differently oriented plane. 
Then we started the dynamics and recorded the behavior of the central spin as a function of time, for different choices 
of Lv and of e. The results are presented in Fig.[5l The central spin can be described by a couple of angles (/), "ip. Here 
we plot the behavior of 4> (similar results are obtained for ?/)) as time passes. In an early stage i < r| the central spin 
remains blocked, since the dynamics starts from the boundary and proceeds towards the interior, as discussed above. 
Then is constant. Later, from t = t-| onwards the central spin begins to rotate unless all the plane is aligned with the 
plane described by the boundary conditions. This is testified by the linear increase of towards a limiting final value. 
The figure shows that the curves for different Ly collapse when time is rescaled as t/Ly. This supports Eq. (|16p . 
Notice also that, by considering different temperatures and e, one obtains data collapse by plotting 0(e/r)^/^. This 
is a consequence of Eq. ^ . Actually di can be written as di = (Se^). Then d is of order \/e{t). Analogously, it can 

be shown that also p ~ yjl{t). Inserting these typical quantities in Eq. ([6|) , and letting cos ~ 1 — , since in the 
late stage spins are rather aligned, one obtains 




In conclusion, the numerical simulation of the evolution of the single domain, confirms our scaling hypothesis. 

At low temperatures, as discussed above, is small and hence the new generated 7 phase is only slightly different 
from the pre-existing (3. Reasoning along the same lines as regarding the simulation of a single domain, the process 
of replacement of an old phase by a new one must be repeated a number UriT) oc T~^l'^ of times in order to 
obtain, in place of the original /3 phase, the phase (say a) of one of the neighboring domains. In conclusion, the 
complete replacement of the phase /3 with a neighboring one (a) requires a time of order At cx nr{T)Ly. When 
this process is completed the typical length of a domain is increased of a quantity ALy(t) oc Lv{t). Therefore 
dLv{t)/dt ~ ALvit)/At oc Lvit)^^T^^^ and hence 

Lv{t) {tT-^f'' (18) 

This prediction can be checked in Fig. H) Here one observes that the curves for Lv{t) relative to different temperature 
quenches collapse (after the pinning) when plotted against tT^^^. The collapse is good for the lower temperatures (for 
T < 10"'*), while it is quite rough at higher temperatures. This is expected since our results are valid in the T —> 
limit. Regarding the power growth law (jlSp . it is satisfactorily confirmed by the data in a certain time window after 
the depinning. For longer times Lv{t) goes abruptly to zero due to the phenomenon of the breakdown of the planes, 
that will be discussed in Sec. IIII Dl As explained in Sec. IIII Dl this phenomenon is delayed lowering T (actually, for a 
relatively high temperature as T > lO""^ it prevents the observation of the law ([18])). 
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FIG. 5: Simulation of a single domain of Lv spins. The angle 4> of the central spin (rescaled by T ' ) is plotted against rescaled 
time t/Ly. 



We consider now the issue of dynamical scaling. In Fig.[6]we plot V{r,t) against x = r/Lv{t) in the range of times 
in which this first dynamical regime occurs. One observes a good data collapse up to r/Lv{t) ~ 1.5 According to 
Eq. this implies that V{r, t) takes a scaling form in this range of r/Lv{t). Since usually scaling is first achieved 
for smaller distances, one could infer that, by pushing the simulations to much longer times, one could observe collapse 
on a larger range of /Lv{t) and conclude that the whole V{r,t) scales. This is surprising, since we have anticipated 
that dynamical scaling is violated in this regime. However this happens because the correlator V{r,t), due to its 
construction, exclusively probes the dynamics of the boundaries of the domains, being blind with respect to the spin 
configuration inside, whose evolution is responsible for the breakdown of dynamical scaling, as we will discuss in 
Sec, nil C 21 One could say that, restricting the attention on the plane boundaries, scaling is obeyed, although globally 
it is not. A similar situation is observed in the d — 1 XY model [2|] where again scaling does not hold (for a different 
reason) but particular correlators, such as C{r,t) or V{r,t) take scaling forms. Clearly this is not a general property 
of every observable, as genuine scaling should imply. Finally, the Porod law (fTO)) is obeyed, signaling that interfaces 
are sharp and that domains remain well defined while coarsening in this whole regime. 



2. Phase- ordering of the spins 

Inside the planes, the spins smoothly rotate generating textures, much like in the one dimensional XY model 
Comparing Fig. [3] and Fig. 2] one understands that the length Lc{t) associated to the spin-spin coherence is much 
smaller then Lv{t), particularly for large times. Then, what really matters for Lc{t) is the evolution of the spins 
in the bulk of the domains, where, since the interfaces are far away, all the additional complications related to their 
presence become irrelevant. As far as the spin spin correlations are considered, therefore, one expects the system to 
behave as a normal conserved vectorial system (but with N = 21): Namely, C{r,t) to obey the scaling form ([5]) and 
Eq. ^ to hold. These features can be checked in Figs. 1317] Regarding the growth law of Lc{t) we obtain a behavior 
in good agreement with what expected, namely Eq. for the higher temperatures, namely for T > 10^*. For the 
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FIG. 6: V{r,t) is plotted against r/Lv{t) for a quench to T = 10 ^ and different times. The dashed hne is the Porod law (|10|l . 



lower temperatures, namely T < 10^^, we measure an effective exponent somewhat larger than 1/4. The case with 
T ~ 10~^ is somehow in between, since the curve initially (after the pinning) grows with an exponent larger than 
1/4 but then the slope is gradually reduced and an exponent in agreement with 1/4 is obtained towards the end of 
the simulation. The behavior of this curve may probably provide an interpretation for what observed for the lower 
temperatures. Namely, the behavior seems to set in after a transient which widens as T is lowered. In the transient a 
slightly larger exponent is observed. Notice that the curves for Lc{t) roughly collapse (we recall that the figure covers 
10 decades in T) when plotted against tT^^^, for the same reason of Lyit). Similarly, the collapse seem to improve 
in quality as T — > 0, as expected, while for T > 10~^ the collapse is poor. 

The scaling form ^ can be checked in Fig. [71 Here one observes a good data collapse up to r/Lc{t) ~ 1.5, similarly 
to what observed for V{r,t). Concerning the shape of C{r,t), differently from V{r,t) it does not display the Porods 
tail, as expected since the rotation of the spins is smooth and there are no sharp interfaces. On the other hand, one 
observes (in the inset) the behavior © typical of vectorial systems but, interestingly, with an effective value N = 2 
which is clearly interpreted as due to the fact that spins in this regime lie on planes. 

The results discussed insofar provide a picture of a system with two different ordering mechanisms at work, coars- 
ening of the planes and phase-ordering of the spins. These profoundly different processes coexist in this regime, 
apparently in a rather independent way, without interfering, possibly because they act on different length-scales. To 
each mechanism a particular correlation function is naturally associated, giving rise to two distinct lengths growing 
with different exponents. Due to that, dynamical scaling is violated even if V{r, t) and C(r, t) possibly scale separately 
with respect to Lv{t) and Lc{t). 
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FIG. 7; C{r,t) is plotted against r/Lc{t) for a quench to T = 10 ^ and different times. The dashed hne is the Porod law (|10|l . 
In the inset, 1 — C{r,t) is plotted against r/Lc{t) on a double logarithmic scale. The dot-dashed line is the form Q, with 
iV = 2. 

D. Breakdown of the domains 

As discussed above, although thermal fluctuations become relevant in depinning the system, their effect in the 
previous regimes is basically to produce the coarsening of the planes without dissolving them. This is because, since 
in the early stage the spins are quite misaligned, both d and p are rather large and hence according to Eq. ([61), the 
typical rotation angle 6 is rather small. 

As the dynamics proceeds, however, textures stretch, spins align, and d and p decrease. In view of Eq. at a 
certain time t^{T), values of 9 sufficiently large, namely cosO — 1 of order unity, become available, which are sufficient 
to destroy the structure of the domains. The breakdown of the planes can be nicely detected by directly looking 
at the spin configuration or, more properly, by inspection of V{r^t). In fact, while for t < Ty{T), V{r,t) takes the 
scaling form (|14p. for t > Ty{T), when the domain disappear, it quickly collapses to a rapidly decaying function. In 
the meanwhile, Lv{t) stops growing and abruptly decreases, as shown in Fig.|4l 

r„ (T) can be evaluated by means of Eq. (fT7|) : The condition cos 6* — 1 ^ 1 for the breakdown of the plains is realized 
when 

e(r,) ~ AT, (19) 

where A is a, constant. This is very well confirmed numerically: In Fig. [8l we plot the average energy computed at the 
time when Ly (t) reaches its maximum, against the temperature of the quench. We find that Eq. (I19p is well verified 
with A ~ 2.76. 

In order to estimate Ty{T) we notice that e(t) is entirely associated with the smooth rotation of the spins in the 
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FIG. 8: The quantity e(r„) is plotted against T. The continuous line is the best fit e{Tv) — AT, with A = 2.76. In the inset 
Ti,(T) is plotted against T. The continuous line is the best Power-law fit Tv{T) — KT~^'^ . 

bulk of the domains, because on the interface between domains spins are perfectly aligned. The typical angle between 
two adjacent spins is ~ L'^^{t). Hence the energy Z{t) — \ — coscj) ^^^{t), a fact that we have verified with good 
accuracy in the simulations. This leads to 

T,{T)^KT-^ , x=l, (20) 

where K is a constant. The dependence of Ty{T) on temperature is shown in the inset of Fig. [51 We find a power-law 
behavior but with a value x — 2-3 only in rough agreement with the expected value x = 5/2. This partial agreement 
is probably due to the fact that our results are only valid in the asymptotic limit of small T and large times t. Here, 
for instance, we cannot consider very small temperatures since for T < 10~^ the time of domain breakdown is too 
long for our simulations. However, as already observed, in such large temperature regime the numerical data do not 
scale according to the asymptotic behavior (recall the discussion regarding Figs l3l4l in Sec. IIII C"|) . 

Since the time to approach and leave the frozen state is negligible in the T — > limit, Tv{T) represents also the 
duration of the regime where scaling is violated. Notice that it increases quite rapidly as quenches are made deeper. 

E. Second phase-ordering regime: absence of domains 

We have seen that the breakdown of the domains occurs when the energy e{t) of the system is comparable to AT 
{A ~ 2.76), Eq. (fT9|l . On the other hand, the system equilibrates at the time Teq{T) when e{t) reaches the equilibrium 
value t{Teq) = Eeq{T) ~ T. Hcncc the energy must still be lowered after the breakdown of the planes. Then, phase- 
ordering must continue even after Ty{T). Clearly, once the domains are eliminated, smooth rotations of the order 
parameter remain the only mechanism at work, and one expects the usual coarsening mechanism of vectorial systems 
for Lc{t), characterized by dynamical scaling with z = 4. This can be checked in Fig. [31 One can observe that the 
power law Q continues to be valid, with no apparent modifications, even after the breakdown of the planes, signaled 
by the decrease of Lv{t) (Fig. (g])). 
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F. Equilibration 

Recalling that Ef.q{T) ~ T, using again e{t) ~ Lc{t)^^ one obtains 

Te,(r) ^ AV„(T). (21) 

This results shows that the duration of this regime characterized by dynamical scaling is comparable to that of the 
previous one. Their duration diverges with the same exponent as T ^ 0. 

In the simulations, after Teg(T) the system is observed to enter the equilibrium stationary state. Computing the 
behavior of some equilibrium quantities, as, for instance, Eeq{T) or ^(T), we found the exact equilibrium results of 
Sec. [IT] with great accuracy. This confirms the correctness and the efficiency of the heat bath transition rates ID 

IV. SUMMARY AND CONCLUSIONS 

In this paper we have studied the kinetics of the one-dimensional Heisenberg model with conserved order parameter. 
The distinguishing feature of this model is the presence of defects in the form of couples of parallel spins separating 
coplanar regions. These arc quite unusual and somewhat counterintuitive defects, since normally one associates the 
notion of defect to regions where the order parameter varies quite abruptly, while here spins are perfectly aligned on 
the defect. Their nature, however, is clearly manifest in the representation of the Vi, where they qualify as unstable 
(but long living in deep quenches), point like defects. Their presence makes the kinetics similar in some respect to that 
of a scalar order parameter, because defects play the role of interfaces in the i7i-representation. In particular, since 
the removal of defects can only be achieved by activated moves, in a low temperature quench the system initially pins; 
the later thermally-activated evolution is characterized by coarsening of the domains with the scaZar-like exponent 
= 3. The vectorial nature of the system makes itself manifest particularly in the smooth rotations of the spins inside 
the domains, producing alignment over a typical length growing with an exponent z = 4 characteristic of vectorial 
systems. This interplay between two different ordering mechanisms continues up to a time Ty{T) (Eq. (I20p ). which 
represents the lifetime of unstable defects. After, defects are removed by thermal fluctuations and a second phase- 
ordering regime sets in, characterized only by smooths variations of the spins, where dynamical scaling is obeyed. 
It is worth mentioning, that the duration of the second phase-ordering regime, without defects, is comparable with 
respect to that of the first one in the low temperature limit. 

These features are unusual and unexpected in non-disordered phase-ordering systems. A natural question, therefore, 
regards their generality, namely if one could expect a similar behavior in other systems. The peculiar dynamics found 
in the Heisenberg chain is obviously related to the conserved character of the kinetics. Therefore we do not expect to 
find something similar in systems without the conservation law, because in that case parallel spins can be singularly 
updated and, in doing so, the defect is removed. In order to check this we have performed simulations of the system 
subjected to the same Hamiltonian but with a dynamics which does not conserve the magnetization. As expected, in 
this case we did not find the unusual features observed with COP, such as scaling violations. The same is found by 
considering a dynamics where conservation is imposed only globally (by exchanging two spins without the constraint 
of neighborhood), as expected since it is known [1„ J_8J that the global conservation law is irrelevant. Restricting to 
systems where conservation is realized locally, as far as we can see there is not a reason preventing the formation of 
similar defects for TV > 3 or, perhaps, even with d > 1. The other ingredient which turns out to be fundamental in the 
model considered insofar is the ability of the system to squeeze all the n spins between two defects into a plane, during 
the regime preceding the pinning. This reduction of the effective internal dimensionality of the order parameter from 
TV to TV — 1 can be achieved because, since n is finite, these spin can be projected on the most energetically favorable 
plane in a finite time (the time over which the system pins at T = 0). On lattices, the feature of n being finite is 
related to the one-dimensional nature of the system. In fact, in d > 1 the typical configuration of the system is a 
bi-continuous percolating structure pL| and, even assuming that a relevant number of couples of parallel spins may be 
formed in a early stage, their geometry should not enclose a domain with a finite number n of spins. Actually, we have 
run some simulation of the Heisenberg model in d = 2 and we have not find the peculiarities of the one-dimensional 
case. On inhomogeneous systems, n is expected to be finite on finitely ramified structures. Therefore, a similar 
behavior could be observed on comb lattices, t-fractals and other finitely ramified networks, where phase ordering for 
discrete models display 1-dimensional features [l^, [l^l . 

On the basis of these reasoning, we also infer that a behavior similar to that of Heisenberg chain with COP could 
be expected for a generic 0{N) model (with iV > 3) with COP in d = 1 or, possibly, on finitely ramified networks. 
Let us mention that our simulations of the 0(4) model quenched to T = show pinning in states characterized by 
defects similar to those of the Heisenberg model. A rather complete analysis of the kinetics for generic iV, similar 
to that presented in this Article, is geometrically rather complicated and beyond the scope of this paper, but may 
represent an interesting issue for further research. 
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